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Microalgae  are  considered  as  a  promising  biofuel  resource  to  fix  carbon  dioxide  (C02)  in  the  future. 
Modelling  microalgae  growth  is  an  effective  method  of  studying  the  performance  of  microalgae  growth 
manifested  by  the  parameters  including  light  distribution,  pigment  dynamics,  nitrogen  uptake,  growth 
rate,  respiration  rate,  temperature  dependence,  and  depth  dependence,  and  helping  to  control  the 
culture  of  microalgae  in  artificial  bioreactors.  In  this  paper,  the  models  of  microalgae  growth  in  nitrogen- 
limited  and  light-limited  culture  system  for  estimating  biomass  productivity  are  investigated  by 
comparing  different  expressions  and  coefficients  used  in  several  main  models.  The  results  show  that 
there  are  some  differences  among  the  results  of  numerical  modelling  due  to  different  expressions  and 
coefficients.  This  paper  will  lay  a  foundation  for  the  selective  use  of  microalgae  growth  models  for 
researchers. 
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1.  Introduction 

2.1.  Characteristics  of  micro  algae 

With  the  continuous  reduction  of  the  fossil  fuels,  and  the 
acceleration  of  global  greenhouse  effect  mainly  due  to  carbon 
dioxide  (C02)  emissions,  energy  supply  may  be  in  trouble  in  the 
near  future.  It  is  urgent  and  significant  to  find  reliable  clean  energy 
resources  alternative  to  fossil  fuels  [1  .  Microalgae  may  become 
one  of  the  most  promising  new  resources  to  supply  energy  and 
mitigate  C02  in  the  future  [2,3].  Technologies  for  producing 
microalgae  and  using  microalgae  for  biodiesel  have  been  known 
for  more  than  50  years.  Up  to  now,  many  researchers  have  done 
plenty  of  work  on  planktons  [4-7  .  Compared  to  the  first  genera¬ 
tion  biofuels,  microalgae  have  several  advantages  in  sustainability, 
economics  and  environment.  Microalgae  not  only  have  higher 
productivity,  but  also  can  be  fed  in  saline/brackish  water/coastal 
seawater  on  non-arable  and  deserted  land.  However,  there  are 
some  shortcomings  of  microalgae  for  use  as  biofuels.  For  example, 
the  dry  solids  content  of  the  microalgae  and  water  mixture  can  be 
as  low  as  0.05%  dry  solids  content.  The  energy  balance  associated 
with  segregating  the  water  from  the  microalgae  dry  solids  to  allow 
a  bio-esterification  process  to  be  undertaken  is  very  high  8]. 
Literature  [9]  also  indicated  that  some  unstable  biodiesel  with 
many  polyunsaturates  will  be  derived  in  the  process  of  producing 
biodiesels  from  algae  oil. 


In  the  above  models,  various  expressions  and  coefficients  have 
been  used.  Investigation  on  the  differences  between  them  is  needed 
so  that  it  is  easier  for  later  researchers  to  select  more  proper 
expressions  and  coefficients  to  model,  predict  and  control  the 
microalgae  growth.  By  now,  little  work  on  the  comparison  of 
different  microalgae  sub-models  (that  is,  the  expressions  or  coeffi¬ 
cients,  for  example,  the  sub-model  of  nitrogen  uptake)  has  been 
reported.  In  this  paper,  based  on  the  valid  Bernard's  dynamic  model 
[20]  acting  as  the  basic  model  (not  the  standard  model)  made  up  of 
several  sub-models,  various  expressions  and  coefficients  in  several 
main  models  [17—19]  will  be  compared  by  comparing  the  basic 
model  and  a  new  model.  The  new  model  is  produced  when  one 
sub-model  of  the  basic  model  is  replaced  by  the  corresponding 
sub-model  used  in  another  model.  In  the  second  and  the  third 
parts,  the  influence  of  other  aspects  such  as  the  depth  of  culture  and 
the  ambient  temperature  is  also  discussed. 

2.  Model  description 

Many  models  have  been  proposed  up  to  now,  which  used 
various  expressions  to  calculate  the  parameters  e.g.  light  distribu¬ 
tion,  pigment  dynamics,  nitrogen  uptake,  growth  rate,  respiration 
rate,  and  temperature  dependence.  In  this  part,  we  will  describe 
Bernard's  model,  and  the  various  expressions  of  parameters  for 
various  models  [17-20]  are  also  discussed. 


2.2.  Microalgae  models 


2.1.  Bernard's  model 


Some  models  have  been  established  by  various  researchers 
to  predict  different  aspects  of  microalgae  growth  process,  such 
as  mass  growth,  light  distribution,  pigment  dynamics,  nitrogen 
uptake,  depth  dependence,  temperature  dependence  and  respira¬ 
tion  rate,  etc.  [10-13  .  The  biological  model  of  microalgae  comes 
from  the  population  model  proposed  by  Malthus  [14  .  The  concept 
of  the  population  model  such  as  growth  rate  is  used  in  later 
microalgae  growth  models.  Droop  [15,16]  proposed  a  dynamic 
model  which  takes  the  dilution  rate  and  influent  inorganic  nitro¬ 
gen  concentration  into  account  to  describe  the  growth  of  micro¬ 
algae.  The  model  is  classical,  but  practical,  which  describes  the 
kinetic  model  of  the  algae  growth  and  the  nitrogen  uptake. 
Recently,  several  new  models  were  proposed  based  on  Droop's 
model.  Geider  et  al.  [17  proposed  a  new  model  which  included 
growth  process,  nitrogen  uptake,  chlorophyll  synthesis,  tempera¬ 
ture  and  respiration  aspects,  but  depth  dependence  was  ignored. 
Quinn  et  al.  [18]  proposed  a  model  for  industrial  scale  systems 
based  on  Geider  et  al.'s  model  [17].  A  dynamic  model  was 
proposed  by  Packer  et  al.  [19]  to  predict  the  growth  and  neutral 
lipid  synthesis  of  green  algae  by  using  a  new  method.  This  method 
considered  the  influences  of  two  factors  (i.e.,  the  photosynthesis, 
and  the  nitrogen  uptake)  on  the  growth  rate,  which  was  defined  as 
a  piecewise  function  of  the  two  factors.  Bernard  [20]  proposed  a 
different  kinetic  model  which  was  verified  by  experimental  data 
[21  .  Different  from  Packer  et  al.'s  model  19],  the  light-limited  and 
nitrogen-limited  factors  will  simultaneously  have  effects  on  the 
growth  rate  in  the  model.  Mairet  et  al.  [22]  further  developed  the 
previous  models  by  dividing  the  biomass  into  three  compartments 
including  the  function  part,  carbohydrates,  and  neutral  lipids.  The 
transformations  and  balances  between  the  three  parts  were  dis¬ 
cussed  in  their  model. 


The  model  [20]  describes  four  variables  in  the  ordinary  diffe¬ 
rential  equations:  s(t),  which  denotes  the  concentration  of  dis¬ 
solved  inorganic  nitrogen,  the  nitrate  or  ammonium  (gNm“3), 
q(t),  which  is  internal  nitrogen  cell  quota  (g  N  (g  C)“  1 ),  x(t)  which 
is  algae  biomass  concentration  (g  C  m-3),  and  /*(t)  which  is  not 
the  real  radiation  but  a  conceptual  variable  denoted  radiation.  The 
unit  of  I*(t)  is  mol  m-2  s-1,  which  means  the  number  of  photons 
absorbed  per  unit  area  per  unit  time.  The  expression  of  calculating 
variable  l*  will  be  discussed  later.  The  four  ordinary  differential 
equations  were  expressed  as: 

s=Ds‘»-<>s*K^-iy-Ds  a) 


q  =  p 


s+K 


x  =  /Z(/0, /*,*,<? ) 


x-Dx-Rx 


i*  =  -p(I0,f,x,q)(  1—^(7—/*) 


where  D  denotes  the  dilution  rate,  p  is  the  maximum  nitrogen 
uptake  rate,  Sin  is  influent  inorganic  nitrogen  concentration,  R  is 
the  inspiration  rate,  I  is  the  average  radiation  along  culture 
volume,  Q0  is  the  minimum  nitrogen  quota,  Q!  is  the  maximum 
nitrogen  quota,  and  f  denotes  the  average  growth  rate,  which  is 
calculated  by 


2k 
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where  J0  represents  the  light  intensity  on  the  bioreactor  surface, 
and  p  is  the  maximum  growth  rate.  The  average  radiation  along 
culture  volume  I  is  given  by 

7  =  ^<l -e~*)  (6) 

where  the  optical  depth,  A,  is  the  product  of  the  depth  of  the 
culture  L  multiplied  by  light  attenuation  rate  £ 

A  =  £L  (7) 

The  light  attenuation  rate  £  can  be  calculated  by 

£  =  aChl  +  bx+c  (8) 

where,  a,  b ,  and  c  are  the  constants,  and  the  chlorophyll  concen¬ 
tration  Chi  can  be  calculated  by 

Chi  =  y(I*)xq  (9) 


The  proportion  of  chlorophyll  concentration  to  nitrogen  con¬ 
centration  y(f)  can  be  given  by 

r(7*)  =  r  (10) 

where  ymax  is  the  maximum  value  of  y(I*). 

The  parameters  A  and  /opfwhich  denote  the  radiation  providing 
maximum  rate  of  photosynthesis  in  Eq.  (5)  are  calculated  by 

A=4I2opt(0co)-K2  (11) 


where  I<si  is 


1KsI 


e\ 


(12) 


(13) 


The  initial  value  of  Chl/x  is  denoted  by  a  parameter 


Chl0 

*o 


(14) 


2.3.  Pigment  dynamics  modelling 

Besides  the  known  a,  b ,  and  c  constants,  the  chlorophyll 
concentration  in  the  reactor  must  be  known  for  calculating  the 
light  attenuation  coefficient.  Bernard  [20]  assumed  that  the 
chlorophyll  concentration  is  proportional  to  the  nitrogen  con¬ 
centration,  xx q,  as  demonstrated  by  Eq.  (9).  The  proportional 
coefficient  is  calculated  with  Eq.  (10).  Quinn  et  al.  18]  considered 
a  =  0  in  Eq.  (8).  It  is  not  necessary  to  calculate  the  chlorophyll 
content  for  calculating  the  light  attenuation  coefficient.  Several 
other  methods  were  also  used  to  calculate  the  chlorophyll 
content. 


(1)  From  Geider  et  al.'s  model 

In  this  model  17],  the  chlorophyll  content  is  calculated  with 
the  following  differential  equation 


1  dChl  _  hch[p 

chl  dt  ~  ec  rchl 


(16) 


where  p  is  the  uptake  rate  of  nitrogen,  0C is  the  ratio  of 
chlorophyll  a  to  phytoplankton  carbon,  rch{  is  the  chlorophyll 
a  degradation  rate  constant,  and  hCM  is  the  chlorophyll  a 
synthesis  regulation  term,  which  is  calculated  by 


=  6 


i  N 

max 


X 

achlOcI0 


(17) 


where  eN  is  the  ratio  of  chlorophyll  a  to  N,  6>j^ax  is  the  maximum 
value  of  0N ,  and  achl  is  the  chlorophyll  a— specific  initial  slope  of 
the  photosynthesis-light  curve. 

(2)  From  Packer  et  al.’s  model 

The  model  19]  used  the  chlorophyll  quota  6C  to  calculate  the 
chlorophyll  content 


(18) 


where  C0  is  the  carbon  subsistence  quota,  and  p  is  the  photo¬ 
synthesis  rate.  Eq.  (18)  didn't  take  the  respiration  into 
account. 


2.4.  Nitrogen  uptake  modelling 


where  Chl0  is  the  initial  value  of  Chl,  and  x0  is  the  initial  value  of  x. 

2.2.  Light  distribution  modelling 

Light  intensity  is  an  important  factor  of  the  algae  growth  model 
[23].  In  nearly  all  the  models,  the  Lambert-Beer  Law  was  used  to 
describe  the  light  distribution  in  the  reactor.  The  law  assumes  that 
the  light  below  the  culture  surface  is  attenuated  at  an  exponen¬ 
tially  decaying  rate  because  the  biomass  and/or  the  chlorophyll 
absorb  a  part  of  the  radiation  from  the  culture  surface 

m  =  I0e~(L  (15) 

Eq.  (6)  can  be  obtained  from  Eq.  (15).  The  exponential  factor  of 
light  attenuation  in  the  culture,  £,  may  be  related  to  the  biomass 
content  and/or  the  chlorophyll  content.  There  are  four  different 
expressions  of  calculating  £:  (1 )  a  ^  0,  b  ^  0,  c  =  0,  that  is,  the  factor 
£  depends  on  both  the  biomass  concentration  and  the  chlorophyll 
concentration,  as  reported  by  [20];  (2)  a  =  b  =  0,  c^0,  that  is, 
parameter  £  =  c  when  the  reactor  depth  L  is  not  too  deep,  as 
reported  by  [17];  (3)  b  ^  0,  a  =  c  =  0,  that  is,  the  light  attenuation 
rate  linearly  correlated  to  biomass  concentration,  as  reported  by 
[18,24];  and  (4)  a  ^  0,  b  =  c  =  0,  that  is,  light  attenuation  coefficient 
£  is  linearly  related  to  chlorophyll  content,  as  reported  by  [11,19]. 


Nitrogen  is  an  important  material  in  microalgae  culture.  The 
nitrogen  source  in  bioreactors  can  be  inorganic  source  such  as 
ammonium  or  nitrate  and  organic  source.  Microalgae  can  synthe¬ 
size  proteins  and  neutral  lipids  by  absorbed  nitrogen.  It  is 
necessary  to  model  the  nitrogen  uptake.  In  Bernard's  model  [20] 
and  Packer  et  al.’s  model  19],  the  nitrogen  uptake  behavior  was 
described  with  Eq.  (2).  Then  a  brief  introduction  to  the  origin  of 
Eq.  (2)  will  be  given. 

The  nitrogen  cell  quota  of  the  culture,  q,  is  equal  to  the  ratio  of 
the  nitrogen  concentration  to  the  carbon  concentration.  This  is 
analogous  to  the  C:N  ratio  which  is  of  crucial  importance  to  the 
digestibility  of  substrates,  q  can  be  expressed  as: 

c  c  v _ y  c  c  y 

<?  =  -=>  <j  =  — Z2^  =  Z~ZS  =  P(S’ q)-KIo,q,Oq  09) 

X  Xz  X  X 

where  s=xxq,  and  p  is  the  growth  rate,  q  changes  due  to  the 
nitrogen  uptake  and  the  microalgae  growth.  The  right  side  of  Eq. 
(19)  can  be  divided  into  two  parts:  the  nitrogen  uptake  part  p(s,q) 
and  the  growth  dilution  part  p(I0,q,^)q 

p(S'q)  =  PsYKs^-^)  (20) 

p(I0,  q,  0  =  Wo,  9  (21) 
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where  p  is  the  maximum  nitrogen  absorption  rate.  Substituting 
Eqs.  (20)  and  (21)  into  Eq.  (19),  Eq.  (2)  is  derived.  Then  several 
methods  to  model  the  nitrogen  uptake  are  shown: 

(1)  From  Geider  et  al.’s  model 

In  this  model  17],  Nitrogen  uptake  is  described  with  a  diffe¬ 
rential  equation. 

q  =  ^-rN  (22) 

q 


uptake  factor  in  Eq.  (27)  linearly  increases  from  0  to  1,  and  that 
in  Eq.  (28)  logarithmically  increases  from  0  to  0.8. 

(3)  From  Packer  et  al.’s  model 

The  growth  rate  is  interestingly  shown  as  a  piecewise  function 
of  time.  The  microalgae  growth  rate  is  determined  by  both 
nitrogen  and  light  factors.  The  growth  rate  is  therefore 
considered  to  be  the  minimum  of  the  growth  rates  calculated 
by  using  the  nitrogen-limited  and  light-limited  relationships 
[19]. 


s  f  Qi -q 
p  Ps+K\  Oi-Oo 


f(T ) 


(23) 


where  rN  is  the  remineralization  rate  constant,  and  f(T ) 
denotes  the  temperature  function. 

(2)  From  Quinn  et  al.’s  model 

In  this  model  [18],  the  temperature  function  f(T )  has  been 
considered.  The  equation  of  nitrogen  uptake  in  Quinn  et  al.’s 
model  is 


where  RN  is  the  respiration  constant  for  nitrogen. 


x  =  px-Dx 


(29) 


i 


H  =  minj  n 


P_1 

q)’c0 J 


(30) 


P  =  QCPm 


(31) 


(qx)2pye/ 
(qx)2  +  (Oo*)2 


(32) 


where  the  value  ofp^,,  is  selected  to  be  7  day-1  in  Eq.  (32),  A  is  the 
optical  cross  section  of  chlorophyll  a,  C0  is  the  carbon  subsistence 
quota,  p  is  3.26  day"  in  Eq.  (30),  and  0  is  the  Quantum  efficiency. 


2.5.  Growth  rate  modelling 


2.6.  Respiration  rate  and  degradation  rate  modelling 


It  is  universally  acknowledged  that  the  growth  rate  is  correlat¬ 
ing  to  the  radiation  intensity.  In  Bernard's  model,  the  nitrogen- 
limited  and  the  light-limited  are  two  multiplication  factors,  as 
shown  by  Eqs.  (5)  and  (21).  Eqs.  (5)  and  (11 )— (14)  are  used  to 
calculate  the  growth  rate  of  microalgae.  There  were  some  methods 
to  calculate  the  microalgae  growth  rate  in  various  models.  Some 
models  considered  the  light  radiation  effects  on  photosynthesis 
rate,  others  (e.g.  Packer  et  al.’s  model  [19]  and  Bernard's  model 
[20])  related  the  growth  rate  to  the  radiation  intensity  directly. 
In  Packer  et  al.’s  model  19],  the  growth  rate  was  determined  by 
both  photosynthesis  rate  and  the  nitrogen  uptake  rate. 

(1)  From  Geider  et  al.’s  model 

Geider  et  al.  [17]  assumed  that  the  growth  rate  is  the  result  of 
subtracting  the  energy  consumption  from  the  photosynthesis 

1  dx 

P  =  xdi  =  p-R^p  ,25) 


There  are  five  sub-models  to  model  the  respiration  or  the 
degradation  of  microalgae.  In  Bernard's  model  [20],  respiration  is 
taken  into  account  only  in  the  growth  modelling.  Quinn  et  al.  [18] 
proposed  a  complex  method  to  distinguishing  two  forms  of  respira¬ 
tion:  respiration  of  growth,  RCl  and  the  reapiration  of  nitrogen  uptake, 
Rn.  The  two  parameters  RN  and  Rc  are  used  in  Eqs.  (24)  and  (25)  to 
represent  the  two  forms  of  respiration. 

In  Geider  et  al.’s  model  [17],  three  respiration  or  degradation 
rates,  chlorophyll  a  degradation  rate  constant,  rch/,  nitrogen 
remineralization  rate  constant,  rN ,  and  Rc  were  used.  The  three 
rates  are  equivalent,  Rc  =  rN  =  rchi  =  0.025  day"1,  which  are  used  in 
Eqs.  (16),  (22)  and  (25),  respectively. 

The  respiration  or  degradation  rate  is  considered  as  a  constant 
in  all  the  above  sub-models.  Mairet  et  al.  [22]  didn't  regard  it  as  a 
constant  respiration  rate,  but  defined  the  respiration  rate  propor¬ 
tion  by  using  a  linear  relationship  to  the  nitrogen  uptake  rate  p 

R  =  r0  +  r^p(s,q)  (33) 


P=Pm 


1 


exp 


achl0cIo 

~P^~ 


(26) 


Pm  =Pr?f(l  \)  (27) 

where  Rc  is  the  maintenance  respiration  rate,  £  is  the  cost  of 
biosynthesis,  pm  is  maximum  value  of  p  at  temperature  T,  and 
p^ef  is  maximum  photosynthesis  rate.  This  model  assumed  that 
50%  of  the  biomass  is  carbon. 

(2)  From  Quinn  et  al.’s  model 

Quinn  et  al.’s  model  18]  is  the  same  as  the  Geider  et  al.’s  one 
[17]  except  the  use  of  different  method  to  calculate  pm 

Pm=P^(l-y)  (28) 


The  nitrogen  uptake  factors  expressions  in  Eqs.  (27)  and  (28) 
are  different.  After  comparing  the  variations  of  the  values  of 
(g-Qo)/(Qi  -Qo)  term  in  Eq.  (27)  and  the  values  of  (1  -Q0/g) 
term  in  Eq.  (28)  with  nitrogen  quota  q  varying  from  0.05  g  N 
(gC)"1  to  0.25  g  N  (g  C)"1,  it  is  found  out  that  the  nitrogen 


where  r0  and  rA  are  the  coefficients. 


2.7.  Depth  dependence  modelling 

In  massive  production,  depth  is  a  vital  parameter  in  bioreactors 
such  as  in  raceway  or  pond.  If  the  bioreactors  are  deeper,  more 
light  radiation  can  be  available  for  use,  but  more  materials  will  be 
wasted.  The  influence  of  the  variable  depth  will  be  discussed  in 
next  part. 


2.8.  Temperature  dependence  modelling 


The  function  related  to  temperature  T  consists  of  two  equations 
which  were  proposed  by  Quinn  et  al.  [25]. 


f(T) 


20  j 

(1  +  0j) 


(34) 


E  E 

0T  =  eR°  T°pt~poT  (35) 

where,  R0  is  the  universal  gas  constant,  Topt  is  the  optimum 
microalgae  growth  temperature,  and  Ea  is  the  activation  energy 
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carboxylation  Rubisco.  Now  the  remaining  problem  is  how  to  add 
the  temperature  function  to  the  ordinary  differential  functions. 
The  temperature  function  can  be  added  to  Bernard's  model  20]  by 
assuming  the  nitrogen  uptake  rate  and  the  average  growth  rate  to 
be  related  to  the  temperature  factor.  The  equations  of  modified 
model  of  Bernard's  are 


s  —  Dsin  p 


s 

s  +  Ks 


f(T)x  -  Ds 


(36) 


q  =  )  (37) 

x  =  p(I0,I*,x,q)(\ -9Af(T)x-Dx-Rf(T)x  (38) 

When  the  environmental  temperature  T  is  equal  to  the  opti¬ 
mum  temperature  Topf,  the  temperature  function /(T)  =  l,  thus  the 
model  will  degenerate  to  Bernard's  model. 

In  Quinn  et  al.’s  model  [18],  the  temperature  factor  is  taken  into 
consideration.  Substituting  temperature  function  into  Eqs.  (20), 
(27)  and  (28),  the  below  equations  are  derived 


3.  Results  and  discussion 

Fig.  1  compares  the  simulations  using  this  model  and  the 
experimental  data  about  Isochrysis  galbana  performed  by  Flynn 
et  al.  [21  .  The  parameters  used  in  the  simulations  with  the  model 
are  presented  in  Table  1.  It  is  assumed  that  no  additional  nutrition 
is  added  into  the  culture  once  the  growth  of  microalgae  began. 
As  shown  in  Fig.  1,  the  model  simulations  have  a  good  agreement 
with  the  experiment  data  [21],  which  shows  the  validation  of  the 
simulations.  This  figure  indicated  that  the  extracellular  nitrogen 
was  exhausted  on  the  12th  day.  At  the  same  time,  the  synthesis  of 
chlorophyll  ceased.  The  growth  process  of  biomass  content  with 
time  could  be  divided  into  four  phases  [21  .  The  first  five  days 
were  the  lag  phase,  during  which  the  cells  increased  at  a  slow  rate. 
Following  the  lag  phase  was  the  exponential  phase  (from  the  5th 
day  to  the  12th  day),  during  which  the  growth  rate  of  cells  was 
0.4  day-1.  The  cells  didn’t  stop  growing  until  the  19th  day,  and 
these  days  (from  the  12th  day  to  the  19th  day)  were  the  post¬ 
exponential  phase.  After  that,  the  culture  entered  the  stationary 
phase,  during  which  the  content  of  biomass,  chlorophyll  and 
intracellular  nitrogen  were  stable. 


<39> 

"--’•'Ai r^)/(T|  (40, 

Pm=P^/(l-y)/C0  (41) 

Based  on  Bernard's  model  [20],  a  modified  Quinn  et  al.’s  model 
[18]  is  established  whose  growth  part  and  temperature  function  are 
substituted  by  Eqs.  (39)-(41)  in  Quinn  et  al.’s  model.  The  model  will 
be  discussed  in  the  next  part. 

Geider  et  al.  [17]  also  took  temperature  effect  into  consideration 
in  nitrogen  uptake,  as  shown  in  Eq.  (23).  Packer  et  al.  [19]  didn’t 
take  the  temperature  factor  into  consideration  in  their  model. 


3.1.  Light  distribution  modelling 

In  order  to  illuminate  the  influence  of  various  light  distribution 
relationships  on  the  simulation  results,  simulations  based  on 
Bernard’s  model  by  using  different  expressions  of  the  light 
attenuation  coefficient  £  are  performed  and  compared. 

The  simulations  show  that,  for  L= 0.05  m,  compared  with  light 
distribution  sub-models  proposed  by  Quinn  et  al.’s  [18]  and  Packer 
et  al.’s  model  19],  the  chlorophyll  values  for  sub-models  from 
Bernard’s  model  [20]  and  Geider  et  al.’s  model  [17]  are  a  bit  lower 
in  algae  biomass  concentration  while  the  inorganic  nitrogen 
concentration,  chlorophyll  concentration  and  intracellular  nitro¬ 
gen  concentration  are  almost  the  same.  Things  are  different  when 
the  reactor  becomes  deeper.  When  L=0.5  m,  the  differences  of  the 
four  curves  corresponding  to  the  four  expressions  of  the  light 
attenuation  coefficient  £  become  obvious  though  their  variation 
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Fig.  1.  Comparison  of  simulation  and  experimental  data  from  [21]. 
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trends  with  time  are  the  same.  Compared  with  light  distribution 
sub-models  proposed  by  Quinn  et  al.’s  [18]  and  Packer  et  al.’s 
model  19],  the  growth  rate  and  nitrogen  uptake  rate  results  for 
sub-models  from  Bernard’s  model  [20]  and  Geider  et  al.’s  model 
[17]  are  lower  while  the  chlorophyll  content  is  even  larger. 

3.2.  Pigment  dynamics  modelling 

Simulations  based  on  Bernard’s  model  [20  by  using  different 
expressions  of  the  chlorophyll  content  are  performed  and  com¬ 
pared.  Fig.  2  shows  that  the  variations  of  four  main  parameter 
values  with  time  calculated  based  on  Bernard’s  model  by  using 
three  different  expressions  of  the  chlorophyll  concentration  Chi 
when  L=0.05  m. 


Table  1 

Parameters  values  used  in  simulation. 


Parameters 

Values 

Units 

P 

1.7 

Day-1 

do 

0.05 

gN(gC)-1 

Qi 

0.25 

gN(gC)-1 

K*si 

1.4 

mol  m  2  s-1 

Ku 

295 

mol  m-2  s-1 

P 

0.073 

g  N  (g  C)_1  day-1 

I<s 

0.0012 

g  N  m-3 

R 

0.0081 

Day-1 

Ymax 

0.68 

g  Chi  (g  N)-1 

kj* 

63 

mol  m-2  s-1 

aChl 

0.46 

m2(gCh/)-1 

/o 

100 

mol  ITT  2  s-1 

C 

2.0 

gC(gN)-1 

r>N 

[/max 

0.3 

g  Chi  (g  N)-1 

tief 

3.0 

Day-1 

C0 

0.61 

gC(gdw)-1 

A 

4.82 

m2g“1  Chi 

< P 

9.84  x  10“2 

g  C  (mol  photons)-1 

Ea 

63 

kj  mol  “ 1 

R0 

8.3145 

J  K-1  mol-1 

T opt 

18 

°C 

S 


As  shown  in  Fig.  2,  the  dissolved  nitrogen  concentration,  and 
the  carbon  and  nitrogen  contents  for  different  pigment  sub¬ 
models  are  almost  the  same  expect  chlorophyll  content.  Fig.  3 
shows  the  four  main  parameter  values  when  L  changes  to  0.5  m. 
The  differences  in  inorganic  nitrogen  concentration,  intracellular 
carbon  concentration  and  intracellular  nitrogen  concentration 
results  of  different  pigment  sub-models  appear  and  the  differences 
between  the  chlorophyll  content  are  still  evident. 


3.3.  Nitrogen  uptake  modelling 

Fig.  4  shows  that  the  variations  of  four  main  parameter  values 
which  is  same  as  Fig.  2  with  time  calculated  based  on  Bernard’s 
model  by  using  three  different  expressions  of  the  nitrogen  uptake 
when  L= 0.05  m.  The  chlorophyll  content,  intracellular  carbon 
concentration  and  intracellular  nitrogen  concentration  result  for 
nitrogen  uptake  sub-models  from  Geider  et  al.’s  model  [17]  is 
smaller  than  the  sub-models  from  Bernard’s  model  [20]  and  Quinn 
et  al.’s  model  [18]  while  the  inorganic  nitrogen  concentration  is 
almost  the  same. 


3.4.  Growth  rate  modelling 

Fig.  5  shows  that  the  variations  of  four  main  parameter  values 
which  is  same  as  Fig.  2  with  time  calculated  based  on  Bernard’s 
model  by  using  four  different  expressions  of  the  growth  rate  when 
L=0.05  m.  Note  that  parameter  p^.ef  is  selected  as  35  day-1  in 
Geider  et  al.’s  model  and  30  day-  in  Quinn  et  al.’s  model. 

As  seen  from  Fig.  5,  it  is  concluded  that  the  four  curves  have  the 
same  trend.  The  growth  rate  for  growth  sub-models  from  Geider 
et  al.’s  model  [17]  and  Quinn  et  al.’s  model  18]  are  drastically 
larger  than  the  sub-models  from  Bernard’s  model  [20]  and  Packer 
et  al.’s  model  [19]  after  the  fourth  day.  But  the  inorganic  nitrogen 
concentration,  chlorophyll  content  and  intracellular  nitrogen  con¬ 
centration  result  for  growth  sub-models  from  four  references 
approximate  to  each  other. 


x 


x*q 


Fig.  2.  Three  different  expressions  to  calculate  the  chlorophyll  content  when  L=0.05  m. 


S.  Yuan  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  34  (2014)  525-535 


531 


3.5.  Respiration  rate  modelling 

Fig.  6  shows  that  the  variations  of  four  main  parameter  values 
calculated  based  on  Bernard’s  model  20]  by  using  five  different 
expressions  of  the  respiration  rate  when  L= 0.05  m.  Based  on 
Bernard’s  model,  the  respiration  parts  in  every  equations  are 
added  in  Eqs.  (2)  and  (3).  The  five  cases  are  analyzed  as  follows: 
(1)  adding  -Rx  into  Eq.  (3),  where  R= 0.0081  day-1;  (2)  assuming 


R  =  0;  (3)  adding  -Rcx  into  Eq.  (3),  where  ftc=4.32xl0-4  day-1; 

(4)  adding  -Rx  in  Eqs.  (2)  and  (3),  where  R=0.0081  day-  ;  and 

(5)  calculating  the  factor  R  with  Eq.  (33). 

As  shown  in  Fig.  6,  the  growth  rate  for  respiration  sub-model 
proposed  by  Mairet  et  al.  [22]  is  a  bit  smaller  than  the  sub-models 
from  the  other  four.  In  general,  the  respiration  factor  is  not  an 
important  aspect  in  the  whole  model,  but  an  accuracy  method  is 
necessary  if  we  want  to  get  a  better  result. 
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Fig.  3.  Three  different  expressions  to  calculate  the  chlorophyll  content  when  L  =  0.5  m. 
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Fig.  4.  Four  main  parameter  values  calculated  based  on  Bernard’s  model  by  using  three  expressions  of  nitrogen  uptake  when  L=0.05  m. 
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3.6.  Depth  dependence  modelling 

In  this  paper,  the  depths  of  bioreactor,  L,  are  0.05  m  and  0.5  m, 
which  represent  the  shallow  bioreactors  and  the  deep  bioreactors, 
respectively.  The  variations  of  the  four  parametric  values  with 
different  depths  are  discussed  and  compared  in  Fig.  7.  As  shown  in 
Fig.  7,  the  results  of  four  main  parametric  values  vary  a  little  with 


the  depth  when  the  depth  L  varies  from  0.03  m  to  0.15  m.  This 
shows  that  it  is  reasonable  to  take  L= 0.05  m  as  a  reference 
example  to  represent  the  shallow  depth  in  this  paper.  Fig.  8  shows 
four  parametric  values  vary  with  time  when  the  bioreactors 
become  deeper.  The  results  of  four  main  values  vary  a  lot  with 
the  depth  when  the  depth  L  varies  from  0.05  m  to  1.00  m.  When 
L=1.00  m,  not  only  the  biomass  concentration  become  lower,  but 


Chi 
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Fig.  5.  Four  main  parameter  values  calculated  based  on  Bernard’s  model  by  using  four  different  expressions  of  growth  rate  when  L=0.05  m. 
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Fig.  6.  Four  main  parameter  values  calculated  based  on  Bernard’s  model  by  using  five  different  expressions  of  the  respiration  rate  when  1=0.05  m.  The  five  expressions  are: 
(1):  adding  -Rx  into  Eq.  (3),  where  K=0.0081  day-1;  (2):  assuming  R  =  0;  (3):  adding  -Rcx  into  Eq.  (3),  where  Kc=4.32xl0-4  day-1;  (4):  adding  - Rx  in  Eqs.  (2)  and  (3), 
where  R= 0.0081  day-1;  (5):  calculating  the  factor  R  with  Eq.  (33). 
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also  more  materials  such  as  nitrate  are  demanded.  Compared 
these  two  figures,  it  is  found  out  that  the  biomass  concentration 
will  decrease  a  lot  when  the  PBR  depth  is  larger  than  0.15  m. 

3.7.  Temperature  dependence  modelling 

Fig.  9  shows  that  the  variations  of  four  main  parametric  values 
calculated  based  on  the  modified  Bernard’s  model  [20]  by  using 


five  different  temperature  rates  when  L= 0.05  m.  As  shown  in 
Fig.  9,  based  on  the  assumption  of  the  optimal  temperature  of 
18  °C  for  Isochrysis  galbana,  the  modified  Bernard’s  model  is 
degenerated  to  Bernard’s  model  without  considering  the  tempera¬ 
ture  factor.  Although  there  are  same  trends  for  all  the  five  curves  at 
various  temperatures  (T=  8, 16, 18,  20,  and  28  °C),  microalgae  grow 
slower  when  the  environment  temperature  is  deviated  from  the 
optical  temperature  of  the  algae. 


Chi 


Fig.  7.  Four  main  parameter  values  calculated  based  on  Bernard’s  model  by  using 


different  expressions  at  five  depths  (L=0.03  m,  0.05  m,  0.08  m,  0.10  m,  0.15  m). 
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Fig.  8.  Four  main  parameter  values  calculated  based  on  Bernard’s  model  by  using  four  different  expressions  at  five  depths  (L=0.05  m,  0.10  m,  0.30  m,  0.50  m,  and  1.00  m). 
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Fig.  9.  Four  main  parameteric  values  calculated  based  on  modified  Bernard’s  model  at  five  different  temperatures. 


4.  Conclusions 

Microalgae  acted  as  the  second  generation  biofuels  could  be 
produced  in  quantity  with  culture  technologies  in  artificial  bior¬ 
eactors.  Modelling  microalgae  growth  is  an  effective  method  to 
predict  the  growth  process,  and  then  estimate  the  production  and 
help  to  improve  the  productivity  of  microalgae  under  a  given 
culturing  condition  at  a  given  site.  Recently,  some  models  have 
been  proposed  with  various  expressions  and  coefficients  of  para¬ 
meters.  In  this  paper,  the  differences  due  to  use  of  the  various 
expressions  and  coefficients  of  parameters  in  the  models  have 
been  compared.  The  conclusions  are  drawn  (based  on  the  case 
L= 0.05  m)  as  follows: 

(1)  Compared  with  light  distribution  sub-models  proposed  by 
Quinn  et  al.’s  [18]  and  Packer  et  al.’s  model  19],  the  chlor¬ 
ophyll  values  for  sub-models  from  Bernard’s  model  [20]  and 
Geider  et  al.’s  model  [17]  are  a  bit  lower  in  algae  biomass 
concentration.  The  inorganic  nitrogen  concentration,  intracel¬ 
lular  carbon  concentration  and  intracellular  nitrogen  concen¬ 
tration  results  of  different  pigment  sub-models  are  almost  the 
same  expect  chlorophyll  content. 

(2)  The  chlorophyll  concentration,  intracellular  carbon  concentra¬ 
tion  and  intracellular  nitrogen  concentration  result  for  nitro¬ 
gen  uptake  sub-models  from  Geider  et  al.’s  model  [17]  is 
smaller  than  the  sub-models  from  Bernard’s  model  [20]  and 
Quinn  et  al.’s  model  [18  while  the  inorganic  nitrogen  con¬ 
centration  is  almost  the  same. 

(3)  The  growth  rate  for  growth  sub-models  from  Geider  et  al.’s 
model  [17]  and  Quinn  et  al.’s  model  [18]  are  drastically  larger 
than  the  sub-models  from  Bernard’s  model  [20]  and  Packer  et 
al.’s  model  19]  after  the  fourth  day.  The  growth  rate  for 
respiration  sub-model  proposed  by  Mairet  et  al.  [22]  is  a  bit 
smaller  than  the  sub-models  from  the  other  four. 

(4)  It  is  concluded  that  the  biomass  concentration  will  decrease  a 
lot  when  the  PBR  depth  is  larger  than  0.15  m.  Microalgae  grow 
slower  when  the  environment  temperature  is  deviated  from 
the  optimum  temperature  of  the  algae. 


This  paper  will  guide  the  researchers  to  select  the  most 
appropriate  models  or  sub-models  or  build  new  models  in  order 
to  model  the  microalgae  growth  process. 

Nomenclature 

A  optical  cross  section  of  chlorophyll  a  (m2 3 4  g- 1  Chi ) 

C0  carbon  subsistence  quota  (g  C  (g  dw)-1) 

Chi  chlorophyll  concentration  (g  Chi  (g  dw)- l) 

D  dilution  rate  (day-1) 

Ea  activation  energy  carboxylation  Rubisco  (J  mol  “  1 ) 

f(T)  temperature  function  (Dimensionless) 

F  nitrogen  restriction  factor  (Dimensionless) 

hchi  chlorophyll  a  synthesis  regulation  term  (Dimensionless) 

70  light  intensity  at  reactor  surface  (mol  m-2  s-1) 

Iopt  radiation  providing  maximal  rate  of  photosynthesis 

(mol  m-2  s-1) 

7*(t)  conceptual  radiation  (mol  m-2  s-1) 

7  average  radiation  along  culture  volume  (mol  m-2  s-  ) 

L  depth  of  culture  (m) 

p  photosynthesis  rate  (day-1) 

pm  maximum  value  of  p  at  temperature  T  (day-1) 

p^ef  maximum  photosynthesis  rate  (day-1) 

q(t )  internal  nitrogen  cell  quota  (g  N  (g  C)_  [) 

Q0  minimum  nitrogen  quota  (g  N  (g  C)-1) 

Q!  maximum  nitrogen  quota  (g  N  (g  C)-1) 

R  respiration  rate  (day-1) 

R0  universal  gas  constant  (J  K-  mol-  l) 

Rc  maintenance  respiration  rate  (day-1) 

rchi  chlorophyll  a  degradation  rate  constant  (day-1) 

rN  nitrogen  remineralization  rate  constant  (day  l) 

Rn  respiration  constant  for  nitrogen  (day-1) 

s(t )  concentration  of  dissolved  inorganic  nitrogen,  (nitrate  or 

ammonium)  (g  N  m-3) 

Sin  influent  inorganic  nitrogen  concentration  (g  N  m-3) 

Topt  optimum  growth  temperature  (°C) 

x(t)  biomass  concentration  (gCm-3) 


S.  Yuan  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  34  (2014)  525-535 


535 


achl  chlorophyll  a-specific  initial  slope  of  photosynthesis- 
light  curve  (m2  (g  Ch/)  —  1 ) 

/(/*)  proportional  coefficient  of  Chlorophyll  concentration: 

nitrogen  concentration  (g  Chl(g  N)- l) 
ymax  maximum  value  of  y(I*)  (g  Chi  (g  N)-1) 

C  cost  of  biosynthesis  (g  C  (g  N) ~ 1 ) 

initial  value  of  6C  (g  Chi  (g  C)-1) 

6C  chlorophyll  a:  phytoplankton  carbon  ratio  (g  Chi  (g  C)-1) 

dN 

maximum  chlorophyll  a:  N  (g  Chi  (g  N)  !) 

6^ax  maximum  value  of  0N  (g  Chi  (g  N)-1) 

A  optical  depth  (Dimensionless) 

ju  growth  rate  (day-1) 

p  average  growth  rate  (day-1) 

p  maximal  growth  rate  (day-1) 

£  light  attenuation  rate  (m-1) 

p  uptake  rate  of  nitrogen  (g  N  (g  C)-1  day-1) 

p  maximum  nitrogen  absorption  rate  (g  N  (g  C)-1  day-1) 

O  Quantum  efficiency  (g  C  (mol  photons)-1)) 
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